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ABSTRACT 

CO measurements of z ~ 1-4 galaxies have found that their baryonic gas fractions are 
significantly higher than galaxies at z=0, with values ranging from 20-80%. Here, we sug- 
gest that the gas fractions inferred from observations of star-forming galaxies at high-z are 
overestimated, owing to the adoption of locally-calibrated CO-H 2 conversion factors (arjo)- 
Evidence from both observations and numerical models suggest that aco varies smoothly 
with the physical properties of galaxies, and that aco can be parameterised simply as a func- 
tion of both gas phase metallicity and observed CO surface brightness. When applying this 
functional form, we find / gas w 10 - 40% in galaxies with M* = 10 10 - 10 12 Mq. More- 
over, the scatter in the observed / gas --W* relation is lowered by a factor of two. The lower 
inferred gas fractions arise physically because the interstellar media of high-z galaxies have 
higher velocity dispersions and gas temperatures than their local counterparts, which results 
in an «co that is lower than the z=0 value for both quiescent discs and starbursts. We further 
compare these gas fractions to those predicted by cosmological galaxy formation models. We 
show that while the canonically inferred gas fractions from observations are a factor of 2-3 
larger at a given stellar mass than predicted by models, our rederived aco values for 2=1-4 
galaxies results in revised gas fractions that agree significantly better with the simulations. 

Key words: galaxies:formation-galaxies:high-redshift-galaxies:starburst-galaxies:ISM- 
ISM:molecules 



is consistent with both observational and theoretical results that 
suggest that even gas rich disc galaxies at z ~ 2 are able to form 
stars rapidly enough that they are comparable to the most extreme 
merger-driven starburst events in the local Universe toaddi et al.l 
120051 ; Irlopkins et al]201ol : lDave et aljioToh . 



However, there is a tension between the inferred gas frac- 
tions of high-z galaxies and galaxy formation models. Hydrody- 
namic cosmological simulations typically account for the simulta- 
neous growth of galaxies via acc retion of gas from the intergalactic 
medium (e.g. iKeres et al]|2003l) as well as the consumption of gas 
by star formation. A generic feature of these simulations is that the 
star formation rate of "main sequence galaxies" (galaxies not un- 
dergoing a starburst event) is roughly proportional to the accretion 
rate, and that galaxies tend to have weakly declining gas fractions 
as their stellar masses increase. Broadly, at a given stellar mass, 
galaxies in simulations have baryonic gas fractions a factor of 2-3 
less than observed gas fractions. This is seen both in hydrodynami c 
simulations, as well as semi-ana lytic models dLagosetalJl201ll) . 
As an example. lDave et al find very few galaxies in a simu- 

lated ~ 150 Mpc (comoving) volume with stellar mass M* ~ 10 11 
Mq with gas fractions greater than 30%. This is in contrast to ob- 



1 INTRODUCTION 



Recent technological advances in (sub)millimetre-wave telescope 
facilities have allowed for the detection of star-forming H2 gas 
in large numbers of galaxies at high-redshift v ia the proxy 



molec ule CO (hereafter, CO; G reve et al. 120051 : iTacconi et al 



20061 : ICoppinetall I2OO8; iTacconi et alj 120081: iBothwell et al 
20 id : iDannerbauer et alj 120091; IWagg et all 120091: ICarilli et al 



iRiechers 



Genzel et al. 2010; Riechers et al. 
20101: ITacconi et all 12010: iGeach et all 1201 ll; 



Casevetal. 2011 



Riechers et al. 201 lab; Iwkng et al.l [2011 



Solo mon & Vanden Bout 1 20051) for a summary of pre-2005 ref- 



erences). 

A major finding from these studies is that, at a given stellar 
mass, early Universe galaxies tend to be significantly more gas rich 
than their present-day counterparts, with baryonic gas fractions-]] 
(hereafter defined as f Ras = Mh2/(Mh2 + M>)) ranging from 
~ 20 - 80% (e.g. Daddi et al. 2010b; Tacconi et al. 2010). This 
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servations which infer gas fractions in comparable mass galaxies 
up to 80%. 

One potential solution is that the inferred gas masses from 
high-z galaxies are systematically too large. H2 masses are typi- 
cally calculated using the luminosity of the CO (J=l-0) emission 
linqfj, and then converted to an H2 mass via a CO-H2 conversion 
factor accEl- In the literature, aco is typically used bimodally 
with one value for "quiescent/disc mode" star formation, and a 
lower value for "starburst/merger mode". In the Galaxy and Lo- 
cal Group, a co is observed to be relatively constant with an aver- 
age aco« 6 telitz et alfeoollFukui & Kawamurall2010h . In con- 
trast, dynamical mass modeling of local galaxy mergers suggests 
that ago should be lower in these galaxie s by a factor of 2-10 
jDownes & Solomon|[l998l ; lNaravanar]|201lh . Despite an observed 
dispersion in inferred aco values from local mergers, a value of 
aco~ 0.8 is typically uniformly applied to these starbursts. 

At higher redshifts, it is more unclear which of the two bi- 
modal values of aco to use. For example, for a star-forming disc 
galaxy that may be undergoing rapid collapse in ~ kiloparsec scale 
clumps and forming stars at rates > 100 MQyr -1 (i.e. compara- 
ble to local galaxy mergers), is the appropriate aco the locally- 
calibrated "quiescent/disc" value, or the "starburst/merger" value? 
Similarly, should high-redshift Submillimetre galaxies (SMGs), 
which are potentially forming stars up to an order of mag- 
nitude faster than local mergers, utilise the locally-calibrated 
"starburst/merger" value? Typically, observational studies use the 
locally-calibrated quiescent/disc value for high-z discs, and the lo- 
cal starburst/merger value for high-z SMGs. 

Recent observ ati onal eviden c e by iTacconi et al. 



Bolatto et alj d2008l); iLerov et alj (1201 ll); iGenzel et alj 



2008); 



2012); 



Schruba et alj d2012t) and Papadopoul os et al] d2012|). as well as 



theoreti c al work by lOstriker & Shettvl J201 lh: INaravanan et al.l 
d2011bh ; IShettv et alj 1201 lal lbh and iFeldmann et alj d2012l) have 
suggested that perhaps the picture of a bimodal aco is too sim- 
plistic, and that aco may depend on the physical environment of 
the interstellar medium (ISM). This picture was expanded upon by 
iNaravanan et alj d2012l) who developed a functional form for the 
dependence of aco on the CO surface brightness and gas-phase 
metallicity of a gal axy. When apply i ng thi s model to observations 
of high-z galaxies. INaravanan et alj J2012b found that on average, 
high-z disc galaxies have aco values a factor of a few lower than 
present-epoch discs, and high-z SMGs have aco values lower than 
present-day ultraluminous infrared galaxies (ULIRGs), with some 
dispersion. Physically, this means that for a given observed CO lu- 
minosity, the inferred H2 gas mass should be systematically less 
than what one would derive using aco values calibrated to lo- 
cal galaxies. This owes to warmer and higher velocity dispersion 
molecular gas in high-z galaxies which gives rise to more CO in- 
tensity at a given H2 column densi ty. The model form for aco 
presented bv lNaravanan et al] j2012t) finds success in matching lo- 
cal observations of discs and ULIRGs (Naravanan et al. 20lTlj), as 
well as observed CO-H2 conversion factors for low metallicity sys- 
tems. 

Building on these results, in this paper, we reexamine CO 



In fact only a few studies directly observe CO (J=l-0) at high-z. Typi- 
cally, higher rotational states are observed, and then down-converted to the 
ground state via an assumption about the CO excitation. 
3 a CO i s alternatively monikered Xqq, or the X -factor. The two are re- 
lated via X C o (cm- 2 (K-kms- 1 )- 1 )=6.3x 10 19 a C o (M Q pc~ 2 (K- 
km s In this paper, we utilise aco as notation for the CO-H2 con- 

version factor. 
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Figure 1. Histograms of literature values for oco f° r high-z galaxies, as 
well as those re-derived via Equation [T] The dashed lines denote the litera- 
ture values, and solid lines the theoretical crco values. The black lines are 
for high-z discs, and red for inferred high-z mergers (SMGs). On average, 
the theoretical ctco values are lower than the locally-calibrated (traditional 
literature values) for discs and mergers. This owes to higher velocity dis- 
persion and warmer gas in high-z galaxies which drives more CO emission 
per unit H2 gas mass than in local galaxies. 



detections from high-z galaxies utilising the physically motivate d 
functional form for aco presented in INaravanan et al. I d2012h . 
rather than the traditional bimodal form. We compare our results 
to those of cosmological hydrodynamic simulations, and show that 
while the inferred gas fractions derived from the traditional aco 
co nversion factor are much larger than those predicted by models, 
the INaravanan et all d2012l) model form for aco brings these val- 
ues down, and in reasonable agreement with simulations. A princi- 
ple result of the work we will present is that the typical gas fraction 
of a high-z galaxy is typically ~ 10 — 40%, ratherthan ~ 40 — 80% 
as is inferred when utilising traditional aco values. In §[2] we de- 
scribe the literature data utilised here; in §[3] we present our main 
results, and in §|4] we summarise. 



2 METHODOLOGY 
2.1 Literature Data 

We examine CO detections of both inferred high-z disc galax- 
ies as well as Submillimetre Galaxies with masses ranging from 
~ 10 10 - 10 12 M Q in stellar masjj . A large n umber of litera- 
ture SMGs are compiled by Bot hwell et alj d2012l) . and include 16 
new detections presented in that pa per. The compila ti on by Both- 
well e t al. includes de t ection s fromlNeri et alj d2 003); Greve et al 
d2005l); ITacconi et al] d200fj); ICasev et alj d2009l) ; iBothwell et al 



J2010h and lEngel etal] d2010h. Oth er SMGs included in our work 
are compiled in IGenzel etal.ld20ld) . The inferred disc galaxies are 
taken primarily from the compilation of IGenzel et alj d2010h and 



4 The limits on stellar masses is highly depen dent on which literature 
stellar masses for SMGs we use iMichatowski et al. 2009; Hainline et al. 
l201ll). 
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Figure 2. Comparison of / ga s against stellar mass (M») for discs and mergers at high-,z. The blue filled circles are observed high-,z discs; the red open squares 
are observed SMGs (assumed to be mergers), and the purple filled circles are optically faint radio galaxies (which are of unknown physical form). The black 
stars show results from main-sequence galaxies fro m the cosmolog i cal hy drodynamic simulations of lDave et alj i2010l) . with dispers ion noted by the error bars . 
The top panels show the results when utilising the Hainline etal. (2011) stellar masses, and the bottom panels when utilising the [Michalowski etal. (2009) 
masses (which are larger by a factor of a few). The left panels show the f gas -M* relation when utilising the traditional locally-calibrated literature aco 
values, and right panel shows the effect of our model «co (Equation[T}. When using the locally-calibrated aco values, the observed gas fractions are a factor 
of 2-3 larger than theoretical gas fractions at a given stellar mass. Because the theoretical aco values are typically lower than the locally-calibrated values 
(FigureflJ, the inferred gas mass for a given CO luminosity decreases, and brings the gas fractions into better agreement with galaxy formation simulations. 



Daddi etall (l2010b|). These inclu de galaxies from the SINS sam- 
ple (Forster Schreiber et al. 2009), as well as SzA"-selected galaxies 
iDaddi et al.1 fe004l) . Finally, we include op tically-faint radio galax- 
ies (OFRGs) with CO detections from the lCasev et alj d201ll) sam- 
ple. 

In our sample, a large number of the BzK and SINS galax- 
ies have been imaged and found to have rotationally dominated 
gas, consistent with a disc-like morphology. The SMGs are often- 
times assumed to be mergers, though there is some debate over this 
jDaveetal.ll2010l ; iNaravanan et alj 120091, l2010bl : lHavward etaD 
l20ld.l201lU2012f) . The OFRGs are of unknown origin. As we will 
discuss in § [3] the global morphology is irrelevant for our model 
form for aco, and the general results in this paper. 

The observational papers that we draw from had to make a 
number of assumptions. Our philosophy is to simply utilise those 
assumptions in this paper, and not make any adjustments to as- 
sumed numbers. The reason for this is to isolate the effects of ap- 
plying our model aco on the inferred gas fractions. For example, 
as we will di scuss, CO surface brightnesses are required in order 
to employ the Naraya nan et alj d2012h model for aco • When direct 
measurements are reported, we utilise those. Otherwise, we make 
the same size assumption that is made in the paper we draw from. 



Similarly, a number of the detections presented in the aforemen- 
tioned papers utilised millimetre-wave telescopes, meaning that the 
observed transition is of higher-lying CO lines in the rest frame. 
Conversion to the ground state CO (J=l-0) line then occurs via an 
assumption of CO excitation. Again, we simply utilise the conver- 
sion from excited CO lines to CO (J=l-0) as presented in the paper 
we pull the data from. This said, the assumed excitation ladders in 
the literature are all relatively similar. 



It is worth a quick word on the stellar masses of the SMGs in 
our sample. There is an ongoing literature debate regarding the stel- 
lar masses of high-r edshift SMGs . Specifically, for th e same SMGs, 
iMichalowski et alj J200gh and lHainline et all d201 ll) find differing 
stellar masses by up to an order of magnitude (with the Hainline 
masses being lower). Some attempts to understand the o r igin of the 
discrepancy have been reported by iMichalowsk iet al.1 J2012h . In 
this work, we remain agnostic as to which stellar masses are "cor- 
rect", and present our results in terms of both sets of observations 
when relevant. 
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2.2 Revised CO-H 2 Conversion Factors for Observed 
Galaxies 

A s discussed in jj[T ] we ut ilise the functional form of aco derived 
in lNaravanan et al.l to re-calculate the H2 gas masses from 

the high- ,z galaxies in our sample. In this model, the CO-H2 con- 
version factor can be expressed as 



10.7 x (Wco) 

£10.65 



(1) 



where aco has units of Mq pc~ /K-km s _1 , Z' is the gas-phase 
metallicity in units of solar, and (Wco) is the luminosity- weighted 
CO intensity over all GMCs in a galaxy. While (Wco) is a difficult 
quantity to observe, in the limit of uniform distribution of luminos- 
ity from the ISM in a galaxy, this reduces to the L' co /A where A 
is the area observed (L' co /A is the CO surface brightness). If the 
light distribution is actually rather concentrated (and most of the 
area observed is in dim pixels), the true surface brightness of the 
pixels which emit most of the light will increase, and the true aco 
will be even lower than what is calculated by Equation[T] This will 
cau se the gas fractions to de crease even further from what utilising 
the [Narayanan e t al.l d2012l) model for aco derives, thus enhanc- 
ing our results. The CO surface brightness ({Wco}) serves as a 
physical parameterisation for the H2 gas temperature and velocity 
dispersion, both of which affect the velocity-integrated CO line in- 
tensity at a given H2 gas mass. 

The functional form for aco also depends on a gas-phase 
metallicity. Physically, aco varies with the gas-phase metallicity 
due to the growth of CO-dark molecular clouds in low-metallicity 
gas. In this regime, the required dust to protect CO from photodis- 
sociating radiation is not present, but the H2 is abundant enough 
to self-shield for survival. For the galaxies in question, metallicity 
measurements are typically not available. Hence, we assume a so- 
lar metallicity (Z 1 — 1) for all galaxies. Based on the z ~ 2 mass- 
metallicity relation, galaxies of mass M, ~ 10 1 1 Mq typically 
have metallicities of order solar dErbet all [2006). Thus, an as- 
sumption of Z' = 1 is likely reasonable. We test the validity 
of this assumption by examining th e effect of includ ing a stellar 
mass-metallicity relation. Following lErb et all ^2006), we assume 
all galaxies above M» = 10 11 Mq have Z' = 1, and that the 
metallicity evolves as (M*) ' 3 for lower mass galaxies. The gas 
fractions in this test do not deviate by more than 10% compared to 
our assumption that Z = This is because of the weak depen- 
dence of metallicity on stellar mass, the weak dependence of aco 
on metallicity, and the fact that gas fractions depend on stellar mass 
as well as gas mass. 

It is worth noting that the results in this paper are not en- 
tirely dependent on the model for aco given in Equation [T] 
A compilation of obse rvations result in a very similar relation. 
lOstriker&ShettvlbOllh showed that aco and £h2 are related via 
a powerlaw in observed galaxies; a conversion of Eh 2 to Wco via 
the relation Xco= £112/ Wco ( an d a linear relationship between 
Xco and aco) gives an aco -Wco relation that is nearly identi- 
cal to the model form in Equation Q] Similarly, while we assume 
Z' = 1 for the observed galaxies analysed in this work, we note 
that the model power-law relation between aco and Z' is very 
similar t o what has been observed in samples of low-metallicit y 
galaxies teolatto et alj2008l ; lLerov et alj201 lklGenzel et alj2012h . 




5 For optically thick gas. 

6 We note that there is one galaxy which va ries by 



22%. This is the 

lowest mass galaxy in the Hainline et al. 1 201 1) stellar mass determinations. 



Figure 3. Difference between observed gas fractions and cosmological sim- 
ulation gas fractions as a function of stellar mas s. The squares are the re- 
sults when util i sing th ejMichaiowski et al. 1 2009) data, and circles from the 
lHainline et alj 1201 ll) data. The open symbols denote the results when the 
observed gas fractions are calculated utilising traditional Qqo values, and 
the fill ed symbols correspond to values calculated from the Narayanan et al. 
120 121) form for aco - O ur functional form for aco results in lower gas 
fractions and, generally, better agreement between the simulations and ob- 
servations. 



In this sense, the forthcoming results in this paper could be derived 
entirely from empirical observational evidence. 

In Figure Q] we plot histograms of the literature aco val- 
ues used in the observations of the galaxies analysed in this work 
(black), as well as our re-derived aco values based on EquationQ] 
(red). We divide the lines into high-z mergers (where SMGs are as- 
sumed to be mergers in this figure and hereafter) and high-z discs 
so that the reader can see the relative difference between our de- 
rived aco values and the original ones. While the range of aco 
values is similar in both cases, there is significantly more power 
toward low aco values when utilising our model for both high-z 
discs and high-z mergers. 



3 RESULTS AND DISCUSSION 

We first examine the effec t of modifying aco fr om the tradi- 
tional bimodal values to the lNaravanan et al] d2012l) model on the 
/g as -M* relation in high-z galaxies. In Figure [2] we show the 
/gas-Aft relation for all gala xies in our sample ut ilising both the 
Michai owski et alj J2009b and lHainline et al.l d201lb stellar masses. 
The left panels show the relationship for the observed galaxies 
when using the traditional bimodal aco, and the right panels when 
applying Equation Q] In order to compare with galaxy formation 
simulations, we overlay the mean / gas - M* relation from the cos- 
mological hydrodynamic calculations o f lDaveetalJ l l2010l) denoted 
by stars. The error bars in the stars denote the range in possible 
/gas values for the simulated galaxies in a given stellar mass bin. 
The simulated galaxies mostly represent "main-sequence" galaxies 
which are not typically undergoing a starburst event. 

When examining the left panels in Figure [2] it is evident that 
the observed galaxies all have substantially higher gas fractions at a 
given stellar mass than the simulations. While lower-mass systems 
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Figure 4. Comparison of / gas against redshift. Filled circles and open squares are same as Figure [2] Green triangles are median and range of data binned 
in bins of Az = 0.5, and are presented to highlight the tre nds of f^ as with redshift. Solid line and dashed line denotes predicted evolution for galaxies in 
(z=0) haloes of M = 10 13 and 1O 14 M0, respectively (Dave et al. 2011). The stars show the mean value and lcr scatter for all simulated galaxies with 
10.5 < log(M«) < 11.5. Both observations and model predict that the mean gas fractions of galaxies should rise with redshift, though utilising our model 
a CO (right panels) versus the traditional literature values (left panels) brings the normalisation of the observed and theoretical / gas -z relation in better 
agreement. 



can be biased to somewhat higher gas fractions because they are 
(in part) selected by far-infrared luminosity, we have tried apply- 
ing similar cuts to simulations and find that this cannot explain the 
difference. In contrast, when applying a CO-H2 conversion factor 
which varies smoothly with physical environment, the inferred H2 
gas masses from the CO line measurements drop and come into bet- 
ter agreement with the simulations. Depending on the stellar mass 
adopted for the SMGs, the gas fractions can drop by up to a factor 
of 3 for a given galaxy. 



The reason for t he drop in gas fraction when using the 
iNaravanan et alj i2012h model for aco versus the traditional bi- 
modal value is due to the typical environments of high- z galaxies. 
The gas fractions of high-z discs are typically large enough that, ab- 
sent substantial internal feedback , large ~kpc-scale clumps of gas 
become unstable and fragm ent iSpringel et al.l200 5 ; Ceverin o et al.l 
l20ld : lHopkinsetal 1 l201 ll ). These clumps can have large internal 
velocity dispersions (~ 50 — 100 kms -1 ), and warm gas tem- 
peratures owing to hig h star formation rates (> 100 MQyr -1 ; 
INaravanan et alj|201 lal) . High velocity dispersions and warm gas 
causes increased CO line luminosity for a given H2 gas mass, 
and reduces aco (Narayanan et al. 201 lb). Because of this, in our 



model, high-z disc galaxies tend to have lowe£] aco values than 
the traditional present-epoch "quiescent/disc" value (though larger 
than the traditional present-epoch "starburst/merger" value; Fig- 
ure QJ. The mean de rived aco for high -z d iscs is 2.5, approxi- 
mately half that of the lDaddi et all feOlObl) and lMagdis eTaUfeOllh 
measurements of high- z BzK galaxies. 

A similar effect is true for high-z starburst galaxies. Owing to 
extreme star formation rat es (potentially up to a thousand MQyr" 1 ; 
lNaravanan" & Dave 2012), the gas temperatures and velocity dis- 
persions in violent z ~ 2 mergers exceed those of even present-day 
ULIRGs. Hence, the average aco is lower than the average ULIRG 
value today. Our average derived value for th e high-z galaxies in 
our sample is aco w 0.5. lMagdiset"al]b01lh finds an upp er limit 
of the aco of a z=4 SMG of 1, and Tacconi et alj j2008l) finds a 
reasonable fit to their observed SMGs with an aco of unity. 

The combined effect of our modeling is that aco for high- 
z discs is typically lower than that of the traditional z=0 "quies- 
cent/disc" value, and aco for high-z starbursts is lower than that 



7 Increased UV photons produced in high SFR galaxies do have the po- 
tential to photodissociate CO. However, these galaxies tend to have large 
dust to gas ratios. Increased dust columns allow GMCs to reach Av ~ 1 
quickl y, and shield CO from photodissociation throughout the bulk of the 
GMC INarav anan et alj2012l) . 
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of the traditional 2=0 "starburst/merger" value (Figure[T|l. Employ- 
ing our model aco consequently lowers gas masses, and brings gas 
fractions into better agreement with cosmological galaxy formation 
models. This is quantitatively shown in Figure[3] where we plot the 
residuals between the observed data and models for both the tradi - 
tional aco , as well as that derived from the lNaravanan et all fc012h 
functional form. 

The usage of our model form of aco reduces the scatter in the 
/gas-M, relation in observed galaxies by a factor ~ 2 at a given 
M» . To calculate the reduction in scatter, we compare the standard 
deviation in galaxy gas fractions within a limited range of stellar 
masses (M, = [5 x 10 10 , lO n ]M0). Much of the scatter in the 
original / gas -M* relation arises from using the bimodal aco val- 
ues. In contrast, our model form of aco varies smoothly with the 
physical conditions in the ISM in a galaxy, and has no knowledge 
as to whether or not the global morphology of a galaxy is a merger 
or a disc. So, if some high-z disc galaxies actually have physical 
conditions in their ISM comparable to starbursts, then their aco 
values will be lower than the canonical "quiescent/disc" aco (Fig- 
ureQ}. The vice-versa is true for high-z SMGs and OFRGs. When 
accounting for the continuum in physical properties in the ISM of 
high-z galaxies (rather than binning them bimodally), the scatter 
in the observed / gas -M* relation reduces substantially. The cor- 
relation coefficient between the observed gas fractions and mod- 
eled ones increases by ~ 10% (from ~ 0.9 to 0.98 for both the 
Michalowski and Hainli ne masses. 

The usage of the iMichaiowski et alj b009l) masses result 
in observed / gas below the simulations for the highest mass 
galaxies. This could reflect either an overestimate of masses, or 
perhap s physical processes that are neglected in the iDave et ail 
(120 1 lh simulations where most SMGs (the most massive galaxies) 
are quiescent, main-sequence objects. Potential neglected phys- 
ical processes include starbu rs ts which may deple te gas (e.g. 
iNaravanan et"ai]|2009l l2010al lbl: lHavward et al]|201lh . or a stage 
of gas consumption without replenishment (that ultimately ends in 
passive galaxies). 

Utilising our model aco additionally results in better agree- 
ment between the observed cosmic evolution of the gas fraction 
of galaxies with redshift and modeled evolution. In Figure [4] we 
plot the observed gas fractions of the galaxies in our literature sam- 
ple against their redshifts. To help guide the eye, we overplot the 
mean values (with dispersion) in redshift bi ns of 0.5. We show the 
predicted values from the analytic model of IDave et al.1 d201ll) for 
halos of mass (at z=0) 10 13 and 10 14 M by th e solid and dashed 
lines respectively, and the simulated points from lDave et alj (2010) 
by the stars. 

When comparing the mean observed values to the predictions 
from analytic arguments and cosmological simulations, we again 
see that when using the traditional aco calibrated to local val- 
ues, the inferred gas fractions are significan tly higher than the pre- 
dictions from models. When applying the Naray anan et al.1 J2012h 
model for aco, the mean values come into better agreement with 
simulations. The scatter (as measured as the standard deviation in 
gas fractions between z — — 3) decreases by ~ 25% when utilis- 
ing our model form for aco as compared to the traditional value 



It is important to note that it is the normal star-forming galax- 
ies (e.g. the BzK galaxies represented by the filled blue circles) that 
come into better agreement with the simulations. On the other hand, 
SMGs, represented by the open red squares in Figure[4] have lower 
gas fractions than the models predict. This is because SMGs are 
not typical galaxies at high-z, but rather rare massive outliers, and 
therefore not reflected in the predictions for an average galaxy at 
high-z. 

Comparing the results of this paper to other galaxy forma- 
tion models is nontrivial. For example, when comparing to galax- 
ies above Lboi> 10 11 Lq , the gas fractions returned from our 
model are significantly lo wer than tho s e pre dicted by recent semi- 
analytic models (SAMs). lLagos et"al] d201ll) utilised the Durham 
SAM to p redict the H2 content in galaxies over cosmic time. As 
shown bv lBothwelletai] J2012h . these models substantially over 
predict the H2 gas fraction as a function of redshift. However, com- 
paring to galaxies above Lboi> 10 Lq produce s better agree- 
ment, however (C. Lagos, private communication). Poppi ng et al.1 
d2012l) utilised an indirect methodology to derive the H2 content 
in observed galaxies. By inver ting the Schmidt relation, and using 
fhe lBlitz & Rosolowskvl 120061) pressure-based prescription for de- 
riving the H2/HI ratio, these authours found a gas fraction-stellar 
mass relation in good agreement with those derived from CO mea- 
surements. Implicit in this model, however, is an assumption of an 
aco conversion factor in setting the normalisation of the observed 
Schmidt relation. In this sense, the measurement is not entirely in- 
dependent of the methods used in CO-derived gas fractions. 



4 SUMMARY 

Observed baryonic gas fractions from high-redshift galaxies as in- 
ferred from CO measurements are typically higher at a given stellar 
mass or redshift than cosmological galaxy formation models would 
predict. These differences can amount to a factor of 2-3 in gas frac- 
tion. 

We suggest that the observed gas fractions are overestimated 
due to the usage of locally-calibrated CO-H2 conversion factors 
(aco)- If aco scales inversely with the CO surface brightness from 
a galaxy (as both numerical models and empirical observational ev- 
idence suggest), then both high-z disc galaxies and high-z mergers 
will have lower average aco values than their z=0 analogs. This 
means that for a given CO luminosity, there will be less underlying 
H2 gas mass. 

Applying a functional form for aco (Equation []} decreases 
the inferred H2 gas masses by a factor of ~ 2 — 3, and brings 
them into agreement with cosmological galaxy formation models. 
Similarly, the usage of our model aco reduces the scatter in the ob- 
served / gas -M* relation by a comparable amount. Galaxy gas frac- 
tions decrease monotonically with increasing stellar mass, while 
the average gas fraction of galaxies in a given stellar mass range 
increases with redshift. 



8 We note that to properly evaluate the evolution of the gas fraction of 
galaxies with redshift, one would ideally examine the same limited stellar 
mass range at each redshift interval. Given the limited number of CO de- 
tections at high-z, however, this is currently infeasible. An examination of 
the galaxies in Figure [2] shows that the majority of our galaxies reside in 
a stellar mass range of log(Af„) = 10.5 — 11.5. The lack of clean sam- 



ple selection is evident in the marginal increase in the already weak corre- 
lation coefficients: the correlation coefficient increases from 0.08 to 0.14 
for the Hainline masses, and remains roughly constant at ~ 0.15 for the 
Michalowski masses. Forthcoming work will address the cosmic evolution 
of galaxy gas fractions in more detail. 
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